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Abstract. 

We study the dynamics of matter waves in an effectively one-dimensional Bose- 
Einstcin condensate in a double well potential. We consider in particular the case when 
one of the double wells confines excited states. Similarly to the known ground state 
oscillations, the states can tunnel between the wells experiencing the physics known for 
electrons in a Josephson junction, or be self-trapped. As the existence of dark solitons 
in a harmonic trap are continuations of such non-ground state excitations, one can 
view the Josephson-like oscillations as tunnelings of dark solitons. Numerical existence 
and stability analysis based on the full equation is performed, where it is shown that 
such tunneling can be stable. Through a numerical path following method, unstable 
tunneling is also obtained in different parameter regions. A coupled-mode system 
is derived and compared to the numerical observations. Regions of (in)stability of 
Josephson tunneling are discussed and highlighted. Finally, we outline an experimental 
scheme designed to explore such dark soliton dynamics in the laboratory. 
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1. Introduction 

One fundamental physical phenomenon observable on a macroscopic scale is the 
Josephson tunneling ( JT) of electrons between two superconductors connected by a weak 
link, predicted by Josephson in 1962 [lj. It is due to the macroscopic wave functions 
with global phase coherence that have a small spatial overlap. The first observation of 
this effect was reported by Anderson et al. [2]. 

Since the only requirement for the occurrence of JT is a weak coupling, other weakly 
connected macroscopic quantum samples were also expected to admit such tunneling. 
For neutral superfluids, JT has been observed in liquid 3 He [I] and 4 He [5]. In the context 
of Bose- Einstein condensates (BECs) [6j [7J [8j [9j [101 E]) the prediction was made by 
Smerzi et al. [T2| [T3| [H] , followed by the experimental observation where a single [151 HH] 
and an array p2] of short Bose- Josephson junctions (BJJs) were realized. The idea of 
BJJs has also been extended to a long BJJ [HI [19], which mimics long superconducting 
Josephson junctions. Such a junction can be formed between two parallel quasi one- 
dimensional BECs linked by a weak coupling. Atomic Bose- Josephson vortices (BJVs) 
akin to Josephson fluxons in superconducting long Josephson junctions [20] have been 
proposed as well [HJ [19] . Moreover, it was emphasized that a B JV can transform from 
and to a dark soliton, due to the presence of a critical coupling at which the two solitonic 
structures exchange their stability. 

The study of JT in BECs considers the tunneling of the Thomas-Fermi cloud, i.e. 
a continuation of the ground state. The tunneling dynamics has been explained using 
a two- mode approximation [121 EH]- The validity of the approximation has been shown 
in [Ml [35] . To improve the applicability regime of such an approximation, modified 
coupled-mode equations have been presented in, e.g., [301 ED E21 133] . 

It is important to note that in addition to the ground state, nonlinear excitations, 
such as dark matter waves, can also be created in BECs. Dark soliton dynamics in 
BECs with single well potentials has been studied theoretically (see a review [22]) 
and experimentally [23l EH [25l [27]. Interesting phenomena on the collective behavior 
of a quantum degenerate bosonic gas, such as soliton oscillations [211 [251 EE] and 
frequency shifts due to soliton collisions [27] were observed. The evolution of solitons 
is of particular interest as the extent to which their behavior can be described in a 
particle picture is an open question and merits further experimental and theoretical 
investigation. A combination of soliton physics with the dynamics at weak links within 
double well potentials will shed light on the collective behavior of excited Bose-Einstein 
condensates in non-trivial potentials. In this paper, we present an analysis of the 
dynamics of dark matter waves in a double well potential. Static properties of such 
a configuration have been recently studied in [281 HZ] • Here, we show that dark matter 
waves can also experience stable quantum tunneling between the wells. This implies 
that localized excitations in higher dimensions, such as vortices, may also experience 
JT. The (in)stability is obtained using numerical Floquet analysis, which is applied 
for the first time in the study of JT. The numerical calculations are necessary as the 
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stability of the observed tunneling is not immediately obvious. This is especially the 
case because dark solitons are higher-order excited states. The possibility that modes 
with lower energy will be excited is not ruled out by a coupled-mode approximation. 

The present paper is outlined as follows. In Section 2, we discuss the governing 
equation used in the current study. We then solve the equation numerically, where 
we obtain stable and unstable Josephson tunneling through a numerical path following 
method. The stability analysis is performed through calculating the Floquet multipliers 
of the solutions. In Section 3, we derive a coupled-mode approximation describing the 
tunneling dynamics. Good agreement between the numerics and the approximation is 
obtained and shown. In Section 4 we present a possible experimental setup to explore 
the results reported herein. Finally we conclude the work in Section 5. 



2. Josephson tunnelings 

2.1. Mathematical model 

We consider the normalized nonlinear Schrodinger (NLS) equation modelling the BECs 
(see, e.g., [29] for the scaling) 

t^t + ^ xx + s\tfj\ 2 ^-V{x)tfj = 0, (1) 

where ip is the bosonic field, and t and x is the time and position coordinate, respectively. 
The parameter s = ±1 characterizes the attractive and repulsive nonlinear interaction, 
respectively, and V(x) is the external double well potential, which for simplicity is taken 

as 

V= 1 -tf(\x\-a) 2 , (2) 

with the parameters Q and a controlling steepness and position of the two minima. The 
total number of atoms N in the trap is conserved with 

/+oo 
\ifj\ 2 dx. (3) 
-oo 

Throughout the present paper, we set s = —1, i.e. we consider repulsive interactions 
between particles. 

For non- interacting particles (s = 0) in a single well potential (a = 0), the governing 
equation (CQ) can be solved analytically to yield ip n = e~ %Ent $ n {x\ where n satisfies 

0n+i = (-y=x - —d x )(j) n , n = 0,1,2, .. ., (4) 

with 

n „.2 



^0 



e 2 ^2 



and the chemical potential E n is given by 
E n = -\/2{2n + l)a 
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The excitations 0o can be continued to nonzero s, which has been considered in, 
e.g., [361 EZl [381 EHl SQl SU S21 S3l HH S51 S6] (see also 08] for discussions on stationary 
solutions of the NLS equation with a multi-well potential that do not reduce to any of 
the eigenfunctions of the linear Schrodinger problem). The existence and the stability 
analysis of continuations of n in a double- well potential has been presented in [2Tj . 
where it was shown that there is a symmetry breaking of the corresponding solutions, 
i.e. a change of (in)stability from a symmetric to an asymmetric state. One typical 
manifestation of the instability is a periodic transfer of atoms between the wells, i.e. 
Josephson tunneling. 

As most of Josephson tunneling studied in BECs considers the tunneling of the 
Thomas- Fermi cloud, which is a continuation of the ground state solution <ft , here we 
consider the tunneling of dark solitons, which can be viewed as continuations of excited 
states n >o- 

2.2. Numerical periodic solutions 

To look for solutions describing Josephson tunneling, we seek solutions that fulfills the 
relation ip(x,T) = ip(x,0), with T being the period of the Josephson oscillations. Such 
solutions posses double periodicity, i.e. one due to the solitonic nature with a period 
2n/E, where E is the chemical potential (intra- well oscillations) and the other one 
caused by the Josephson effect (inter- well oscillations). Consequently, we can express 
the solutions in terms of a Fourier series multiplied by a factor related to the stationary 
character of dark solitons 



where u = 2ir/T is the Josephson oscillation frequency. These solutions are denoted as 
commensurate if the commensurability condition E = (q/p)T = (2qn)/(ojp) is fulfilled, 
with {q,p} G N. In what follows, we fix p = 1. 

Commensurate solutions are consequently fixed points of the map ip(x, 0) — > ip(x, T) 
and can be found either by using shooting methods in real space or algebraic methods 
in Fourier space. In order to do that, we will transform the problem into a discrete 
one by means of a finite difference discretization with spatial step Ax = 0.2 and apply 
the techniques developed for discrete breathers in Klein-Gordon lattices [50] . If a 
shooting method were used, a time step At = 0.02 would be necessary. As the considered 
oscillations herein have periods about 1500 time units, this method would imply many 
integration steps. In addition to that, the lack of an analytical Jacobian would also 
imply the necessity of the numerical determination of this matrix. These facts suggest 
the suitability of the proposed Fourier space method, which, apart from transforming 
the set of differential equations into an algebraic one, provides an analytical expression 
for the Jacobian. 

Truncating the Fourier series at k m , i.e. the maximum value of \k\, which has been 
chosen to be 9 in most of the calculations due to computational reasons, Eq. ([!]) yields a 



oo 




(5) 



fc=- 



oo 
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set of nonlinear equations with the k-th component of the dynamical equation set given 
by 

F k (x) = (E — uk)z k + d\z k - V(x)z k -s ^ ^2 Z mZ n Z k - m +n = 0.(6) 
We then obtain the following expression for each component of the Jacobian 



dF k (x) 
dzJx') 



{[E-cok- V{x))5{x - x') + d 2 xx ] S k>n 

— s5(x — x ) [z m z k - n+m + z m (z k - m+n + z n+m _ k )\ , (7) 



where we have written z k = z k (x) in both equations. 

Once a periodic solution, say *ff(x,t), is obtained, to study its (linear) orbital 
stability one needs to analyze the time evolution of a small perturbation £ (x, t) to 
*&(x, t). The equation satisfied to leading order by £(x, t) is 

%t + U - s(m 2 Z + ^ 2 D - V(x)£ = 0. (8) 

The stability properties are then determined by the spectrum of the Floquet operator 
J 7 (whose matrix representation is the monodromy) defined as 



Re(£(x,T)) 
Im(£(x,T)) 



(9) 



Re(£(z,0)) 
Im(£M)) 

A Floquet analysis can be performed as long as the solutions are commensurate. 

As the Floquet operator is symplectic, it implies that there is always a pair of 
degenerate monodromy eigenvalues corresponding to the phase and growth modes at 1. 
If the oscillations are stable, all the eigenvalues must lie on the unit circle (see [51] for 
a similar analysis in a discrete setting). In order to get the monodromy with enough 
accuracy, the simulations must be performed using a time step twenty times smaller 
than in the case of the dynamical equations, i.e. At = 0.001 in this case. 

We have calculated commensurate solitons for Q = 0.1 and a = 10 using the 
method described above and analysed the stability of those solutions. Presented in the 
top panels of Figure [U are two periodic solutions that we obtained together with the 
time evolution of a dark soliton in a double well potential. The left and right panel 
respectively corresponds to JT and a transition to macroscopic quantum self-trapping, 
similarly to the dynamics of the ground state oscillations (T2J UJ] . 

In the middle panels of Figure [fl we present the distribution of the Floquet 
multipliers of the two solutions depicted in the top panels in the complex plane. Note 
that we did not obtain a continuum spectrum of the Floquet operator due to the 
discretization of the equations. It is worth noting that as there is a quartet of multipliers 
that do not lie on the unit circle, one can conclude that the solution in the top right 
panel is unstable. We show in the bottom panels of Figure [1] a longer time evolution 
of the solutions in the top panels, where one can see that the solution in the top right 
panel is indeed unstable. The instability we reported here is a clear evidence that 
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Figure 1. (Top) The first few oscillations of the atom density \tp{x,t)\ 2 for dark 
solitons in a double well potential with O = 0.1, a = 10, and (left) uj — 0.00450 and 
(right) oj = 0.00520, which respectively corresponds to N = 0.0340 and N = 0.7677. 
In both cases, the initial conditions are obtained from a numerical continuation with 
q = 47 (see the text). (Middle) Floquet multiplier distributions corresponding to 
solutions in the top left and right panel, respectively. (Bottom) Longer time evolutions 
of the top panels where one can see that the solution in the top right panel is indeed 
unstable. 
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the nonlinearity term in the governing equation (pQ) plays an important role, as all 
the solutions would have been stable otherwise. A typical instability dynamics is a 
repulsive interaction between the dark solitons in different wells so that they start to 
oscillate about the minimum of the wells. This can be clearly observed in the bottom 
right panel of Figure [U 

We have also obtained periodic solutions for various parameter values. In the top 
left panel of Figure [2] we show the dependence of the norm (number of atoms) N of 
tunneling dark solitons when the inter-well oscillation frequency is varied. In the panel, 
several representative values of q are considered. Note that the possible values of q are 
not limited to those shown in the graph. As uj is increased further, there is a critical 
value above which solutions are unstable. Unstable solutions are indicated as dashed line 
in the top left panel. The solutions can also be continued for decreasing frequencies u 
down to a critical value for which the solutions transform into a non-oscillating one (not 
shown here). In the top right panel of Figure [2] we show the dependence of the growth 
rate (the logarithm of the maximum modulus of the Floquet operator eigenvalues) with 
respect to u for q = 47. We also present the growth rate of JT for a fixed u and q and 
varying separation distance between the two wells a in the bottom panels of the same 
figure, i.e. u = 0.0049 and q = 47. For small a, the solutions tend to a non-oscillating 
one with one dark soliton in each well, analogously to what occurs for small u and fixed 
a. 



3. Coupled-mode approximations and their validity 

To describe dark soliton dynamics reported in the previous section, we will readily use 
a two-mode approximation derived in [31]. Following [31] . we write 

V, = VN (b 2 (t)<f> 2 (x) + &s(t)S 3 (aO) > $ 2,s = S+fr^-fr) , (io) 

where $±(x) is a continuation of 02,3 © for nonzero a satisfying 

d m & ± + /3±$± - V(x)& ± + sN$l = 0, (11) 

with the constraint J &j&k dx = 8j,k, hi = +, — • Two examples of $j, which 
corresponds to the norm iV in the Figure [1] are presented in Figure El 

Next, for simplicity we write bj(t) = \bj(t)\e t0i ^ . Equations ([3]) and ([TO]) imply that 
\b 2 (t)\ 2 + \b 3 {t)\ 2 = 1. Defining 

z(t) = \b 2 (t)\ 2 - \b 3 (t)\ 2 , A6(t) = 9 3 (t) - 9 2 (t), (12) 

one can obtain the equations satisfied by z and A# [31] 
dz _dH dA6 _dH 

~dl~~^Ke 1 ~dT~~dz~' ( ' 

where 

H = ^Az 2 - BVT^z 2 cos Afl + -C(l- z 2 ) cos 2Afl, (14) 
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Figure 2. The top left panel presents the dependence of the norm with respect to 
lu for dark solitons with a = 10. Dashed lines indicate unstable solutions. Here, q 
sweeps the values between 40 and 50. The top right panel shows the dependence of 
the growth rate with respect to uj for q — 47. Bottom panels depict the norm and the 
growth mode of tunneling dark solitons with fixed ui = 0.0049 and q = 47 for varying 
a. 



A 
C 



IO7+- - 7++ - 7- 



B 



-27+- + 7++ + 7— 



P--P+ + 
ljk = -sN 



7++ - 7- 



(15) 

(16) 
(17) 

with j, k = +,—. 

We plot the phase-portrait of ( fl~3l) in Figure H] for the two values of N in Figure [Q 
To compare the two-mode approximation with the top panels of Figure [H we calculate 
the variable z from the numerics of the full equation ([1]) as [31] 



f°_J^t)\*dx-N/2 
NS 



S 



$ , $_ dx 



where in the present case S ~ 0.5. As A6 can be calculated immediately, one can 
compare the numerics and the approximation right away. Shown in Figure S] are the 
comparisons, where satisfactory agreement is obtained. As for the instability of the 
solution in the top right panel, that develops at a later time, it is beyond the validity 
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Figure 3. The second and third collective modes of the confining potential V(x) 
(dash-dotted) for (solid) N = 0.034 and (dashed) N = 0.7677, with s = -1. 
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Figure 4. The phase-portraits of (|13|) for the two values of N in Figure Q] i.e. (left) 
N = 0.0340 and (right) N = 0.7677. Thick symbols correspond to the periodic 
solutions shown in Figure [T] 



of any currently available two-mode approximations. One needs a better ansatz for the 
approximations to capture the stability of the periodic solutions. Note that the validity 
issue discussed herein is completely different from that in [31]. In [31], the issue is 
related to the fact that the approximation does not capture the Josephson oscillation 
of the full equation directly from the beginning, which typically occurs when \sN\ ^> 1, 
while in our case \sN\ < 1 and the approximation does capture the existence, but not 
the stability. 

One can observe that the phase portrait in the left panel of Figure H]has two families 
of periodic oscillations, i.e. one centred at A6 = and the other at ±7r. The latter is 
known as 7r-oscillations [13]. The stable solution in the top left panel of Figure [2] with 
q = 50 and the same norm belongs to this family. With this, we conjecture that all 
trajectories in the phase portrait are stable. 
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For the phase portrait in the right panel of Figure HJ one can observe that there are 
two types of solutions, i.e. Josephson oscillations and running states. The latter type is 
also referred to as macroscopically quantum self-trapped states. As we have seen that 
the two-mode approach can provide a good approximation to periodic solutions, one can 
use the ansatz (flOl) and direct simulations of ([I]) to predict whether or not a periodic 
solution is stable. With this, we conjecture that all trajectories in the phase portrait 
are unstable. 

Based on our analysis above, it is likely that all periodic solutions belonging to the 
same phase portrait, i.e. having the same norms, will have the same stability property. 
Nonetheless, one can easily notice in the top left panel of Figure [2] that the critical 
norms, above which the corresponding solutions are unstable, for different values of u 
are not exactly the same. Further analysis whether or not the discrepancy is caused by 
the finite number of Fourier modes k m is to be addressed in future investigations. 

4. Experimental setup 

Dark solitons can be created in a magnetic trap in various ways, including phase- 
imprinting [231 EH EZ] , merging of two condensates [25], and passing a penetrable barrier 
through condensates [52]. In the context of condensate splitting and merging, a method 
of applying oscillating radio- frequency (rf) fields in combination with static magnetic 
fields has been recently proven to facilitate good experimental control over creating, 
tuning, and manipulating double well potentials [531 EH E5]. Without changing the 
static field configuration that provides to very good approximation a parabolic single 
well trap, for example in atom chip based microtraps [56], modifications of frequency 
and amplitude of the rf field allow for raising and lowering a splitting potential barrier. 
The barrier height between the wells is readily controlled in this type of setup. Moreover, 
inhomogeneities of the rf field, introduced by locally producing the field by an additional 
conductor on an atom chip [53], can be exploited to introduce an imbalance or slight 
asymmetry between the two wells. Note that rf field engineering in the context of 
microtraps with typical distances between trapped atoms and field sources on the order 
of 10 — 100 (ttm is straightforward as the wavelength of the oscillating field (in the MHz 
range) is by far sufficiently long to warrant a near-field treatment, with a DC calculation 
yielding accurate results. 

Figure [5] illustrates a possible scheme based on these methods to produce a soliton in 
an originally harmonically trapped degenerate Bose gas and to observe its subsequent 
Josephson oscillations in a double well potential. The protocol starts with a simple 
harmonic magnetic trap, as can be produced on an atom chip by a Z-shaped wire [56J. 
The gas is then split with rf fields into a symmetric double well. The double well is then 
slightly imbalanced so that after a time of typically a few ms, a relative phase difference 
of the two clouds of <ft — 77 will have accumulated. Merging the potentials will now 
result in a solitonic excitation of the combined gas (see, e.g., [57]). Slight controllable 
deviations A0 from it will produce a (slowly) oscillating soliton. After raising the barrier 
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Figure 5. Schematic implementation of a protocol for controlled production of a 
soliton in a double well potential (see text). 

once again, the tunneling dynamics of the soliton can now be studied. The scheme can 
be extended to multiple solitons by splitting the potential in more than two wells, for 
example by applying multiple frequency component rf fields [58], relevant in the study 
of interactions of different numbers of dark solitons in each well. 

5. Conclusion and future work 

We have studied dark soliton dynamics in a double well potential, where it has been 
shown that dark solitons can experience JTs between the wells. A coupled-mode 
approximation has been derived to explain the observations. Numerical stability analysis 
based on the full governing equation has been performed to show that JTs can be 
stable. Through path following methods, unstable solutions were also obtained. An 
experimental scheme designed to explore such dark soliton dynamics in the laboratory 
also has been outlined. 

A natural problem to follow the tunneling reported above is when two different 
numbers of dark solitons are loaded into each of the minima. We present in Figure 
|6] interactions between one and two dark solitons. One can observe that there is an 
interference pattern analogous to the acoustic beating pattern in the interaction of two 
continuous waves with slightly different frequencies. As shown in the left panel indicated 
by the dashed ellipses, there are two levels of modulated patterns; the big ellipse shows 
one tunneling period modulated by the oscillation in the small ellipse. In the right 
panel, we zoom in on the small ellipse to show that a beating pattern also occurs on a 
smaller scale. It can be calculated that the oscillation period in the small ellipse when 
a ^> 1 is approximately T w 2ir/Q, which is in accordance with the numerical result. A 
multi-mode approximation can be obtained as before as briefly discussed in [31]. It is 
then interesting to study the stability of such interactions. Together with the question 
on a better ansatz that can predict the stability of periodic solutions reported herein, 
this is currently under investigation and will be reported elsewhere. 
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Figure 6. (Color online) Interactions between two different excitations when loaded 
into each of the minima of the magnetic trap. Big and small dashed ellipse indicates 
one tunneling period of the excitation in the left and right well, respectively. The 
bottom panel zooms in the small dashed ellipse. The total number of atoms is small, 
similar to that in Figure [TJ 



Acknowledgment s 

We acknowledge fruitful discussions with Andrea Trombettoni and Panos Kevrekidis. 
JC acknowledges financial support from the MICINN project FIS2008-04848. 

References 

[1] B.D. Josephson, Phys.Lett. 1, 251 (1962). 

[2] P.L. Anderson and J.W. Rowell, Phys. Rev. Lett. 10, 230 (1963). 

[3] A. Barone and G. Paterno, Physics and Applications of the Josephson Effect (John Wiley & Sons, 

New York, Singapore, 1982). 
[4] S.V. Pereverzev, A. Loshak, S. Backhaus, J.C. Davis, and R.E. Packard, Nature 388, 449 (1997). 
[5] K. Sukhatme, Y. Mukharsky, T. Chui, and D. Pearson, Nature 411, 280 (2001). 
[6] S.N. Bose, Zeit. Phys. 26, 178 (1924). 

[7] A. Einstein, Sber. Preuss. Akad. Wiss. 22, 261267 (1924); ibid. 1, 314 (1925). 
[8] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E.A. Cornell, Science 269, 
198201 (1995). 

[9] K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn, and W. 
Ketterle, Phys. Rev. Lett. 75, 39693973 (1995). 
[10] C.C. Bradley, C.A. Sackett, and R.G. Hulet, Phys. Rev. Lett. 75, 16871690 (1995). 
[11] F. Dalfovo et al., Rev. Mod. Phys. 71, 463512 (1999). 

[12] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Phys. Rev. Lett. 79, 4950 (1997). 
[13] S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. A 59, 620633 (1999). 
[14] S. Giovanazzi, A. Smerzi, and S. Fantoni, Phys. Rev. Lett. 84, 4521 (2000). 
[15] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M.K. Oberthaler, Phys. Rev. Lett. 
95, 010402 (2005). 

[16] S. Levy, E. Lahoud, I. Shomroni, and J. Steinhauer, Nature 449 579 (2007). 

[17] F.S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A, Trombettoni, A. Smerzi, and 

M. Inguscio, Science 293, 843 (2001). 
[18] V. M. Kaurov and A. B. Kuklov, Phys. Rev. A 71, 011601 (2005). 
[19] V. M. Kaurov and A. B. Kuklov, Phys. Rev. A 73, 013627 (2006). 



Josephson tunneling of dark solitons in a double-well potential 



13 



[20] A. V. Ustinov, Physica D 123, 315 (1998). 

[21] G. Theocharis, P.G. Kevrekidis, D.J. Frantzeskakis, and P. Schmelcher, Phys. Rev. E 74, 056608 
(2006). 

[22] V.V. Konotop, in Emergent Nonlinear Phenomena in Bose-Einstein Condensates Theory and 

Experiment, Eds. P.G. Kevrekidis, D.J. Frantzeskakis, R. Carretero-Gonzalez (Springer Berlin 

Heidelberg, 2008) pp. 65-97. 
[23] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sengstock, A. Sanpera, G. V. Shlyapnikov, and 

M. Lewenstein Phys. Rev. Lett. 83, 5198 (1999). 
[24] C. Becker, S. Stellmer, P. Soltan-Panahi, S. Dorscher, M. Baumert, E-M. Richter, J. Kronjager, 

K. Bongs, and K. Sengstock, Nature Phys. 4, 496 (2008). 
[25] A. Weller et al., Phys. Rev. Lett. 101, 130401 (2008). 
[26] G. Theocharis et al., Phys. Rev. A 81, 063604 (2010) 

[27] S. Stellmer, C. Becker, P. Soltan-Panahi, E.-M. Richter, S. Dorscher, M. Baumert, J. Kronjager, 

K. Bongs, and K. Sengstock Phys. Rev. Lett. 101, 120406 (2008). 
[28] S. Middelkamp, G. Theocharis, P. G. Kevrekidis, D. J. Frantzeskakis, and P. Schmelcher, Phys. 

Rev. A 81, 053618 (2010). 
[29] R Carretero-Gonzalez et al., Nonlinearity 21, R139 (2008). 

[30] E. A. Ostrovskaya, Yu. S. Kivshar, M. Lisak, B. Hall, F. Cattani, and D. Anderson, Phys. Rev. 
A 61, 031601-4 (2000). 

[31] D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 (2006); Phys. Rev. A 74, 039905(E) 
(2006). 

[32] XY. Jia, WD. Li, and J.Q. Liang, Phys. Rev. A 78 023613 (2008). 

[33] B. Julia-Diaz, J. Martorell, M. Mele-Messeguer, and A. Polls. larXiv:1012.1023l 

[34] A. Sacchetti, SIAM J. Math. Anal. 35, 1160-1176 (2003). 

[35] A. Sacchetti, J. Evol. Eq. 4, 345-369 (2004). 

[36] J. A. Vaccaro, O. Steuernagel, and B. Lorimer, ArXiv:cond-mat/0007233vl. 

[37] M.P. Coles, D.E. Pelinovsky, and P.G. Kevrekidis, Nonlinearity 23, 1753-1770 (2010). 

[38] R.J. Dodd, J. Res. Natl. Inst. Stand. Technol. 101, 545 (1996). 

[39] Yu.S. Kivshar, T.J. Alexander, and S.K. Turitsyn, Phys. Lett. A 278, 225 (2001). 

[40] V.I. Yukalov, E.P. Yukalov, and V.S. Bagnato, Phys. Rev. A 56, 4845 (1997). 

[41] C. Trallcro-Gincr, J.C. Drake-Perez, V. Lopez-Richard, J.L. Birman, Physica D 237, 2342 (2008). 

[42] P.G. Kevrekidis, V.V. Konotop, A. Rodrigues and D.J. Frantzeskakis, J. Phys. B: At. Viol. Opt. 

Phys. 38, 1173 (2005). 
[43] G.L. Alfimov and D.A. Zezyulin, Nonlinearity 20, 2075 (2007). 

[44] D.A. Zezyulin, G.L. Alfimov, V.V. Konotop and V.M. Perez-Garci'a, Phys. Rev. A 78 013606 
(2008). 

[45] R. D'Agosta, B. A. Malomed, and C. Presilla, Laser Phys. 12, 37-42 (2002). 
[46] R. D'Agosta, B.A. Malomed, and C. Presilla, Phys. Lett. A 275, 424-434 (2000). 
[47] R. Ichihara, I. Danshita, and T. Nikuni, Phys. Rev. A 78, 063604 (2008). 
[48] R. D'Agosta and C. Presilla, Phys. Rev. A 65 043609 (2002). 
[49] J. L. Marin and S. Aubry. Nonlinearity 9, 1501 (1996). 

[50] J. Cuevas, J.F.R. Archilla and F.R. Romero. J. Phys. A: Math, and Theor. 44, 035102 (2011). 

[51] T. R. O. Melvin, A. R. Champneys, P. G. Kevrekidis and J. Cuevas. Physica D 237, 551 (2008). 

[52] P. Engels and C. Atherton, Phys. Rev. Lett. 99, 160405 (2007). 

[53] T. Schumm et al, Nature Physics 1, 57 (2005). 

[54] I. Lesanovsky, Phys. Rev. A 73, 033619 (2006). 

[55] S. Hofferberth et al., Nature Physics, 2, 710 (2006). 

[56] R. Folman et al., Adv. At. Mol. Opt. Phys. 48, 263 (2002). 

[57] A. Negretti and C. Henkel, J. Phys. B: At. Mol. Opt. Phys. 37, L385 (2004). 

[58] I. Lesanovsky and W. von Klitzing, Phys. Rev. Lett. 99, 083001 (2007). 




3000 
2500 
2000 
1500 
1000 
500 




-20 20 



